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We study the critical behavior of the single-site entanglement entropy S at the Mott 
metal-insulator transition in infinite-dimensional Hubbard model. For this model, the 
entanglement between a single site and rest of the lattice can be evaluated exactly, 
using the dynamical mean- field theory (DMFT). Both the numerical solution using 
exact diagonalization and the analytical one using two-site DMFT gives S — S c oc 
alog 2 [(1/2 — D c )/D c ] (U — U c ), with D c the double occupancy at U c and a < being 
different on two sides of the transition. 
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1. Introduction 

The concept of quantum entanglement plays a key role in the field of quantum 
information manipulation and processing^. It describes the inseparability between 
parts of a given system in a given state. In the past decade, the close relation 
between quantum entanglement and the quantum correlation in many-body systems 
are discusseoP, especially in the context of condensed matter physics^, cold atomJ^l, 
and quantum chemistry^, etc. 

One of the most interesting ideas is to employ the entangleme nt entropy as an 
indicator of quantum phase transiti on (_QPT j~j n boi li .-pin systema ^^^^^^^^^^^ l 
and the interacting fermion system a ' ' ' ' ' ' '. For those quantum phase 
transitions that cannot be described by Landau's symmetry breaking paradigm, 
there is no well defined local order parameters. In such cases, the entanglement 
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entropy may be a really useful concept for characterizing the quantum phase. 

One of such QPT is the Mott metal-insulator transition, where a many-body 
system transits from a metallic state into an insulator, when the interaction strength 
between particles exceeds a critical value^H Experimentally, Mott transition (MT) 

|O0| 

has been widely studied both in strongly correlated electron systems such as V2O3 
and in cold atom systems^ 1 . In the most strict sense of the Mott transition, no 
symmetry breaking occurs at the transition and hence it belongs to the type outside 
the Landau's paradigm. Since 1960's, Mott transition has been one of the key 
issues in condensed matter physics. Theoretically, intensive studies based on the 
Hubbard-type models have been done in the past decades^l. The study of MT in 
terms of the quantum entanglement, e specially usi ng the entanglement entropy as 
a measure, appears for one- dimensional™™] two-dimensionai^, and infinite 
dimensional^ Hubbard-like models. 

For fermionic lattice models, a connection has been established rigorously be- 
tween the singularity of single-site entanglement entropy S and the order of QPT, 
under certain conditions^! This connection states that the discontinuity in the (k 
- l)-th order derivative of S gives a k-th order QPT. For one dimensional Hubbard 
model, 5* reaches a maximum at U c — where MT occur due to equal popu- 
lation of all the local bases at the transition poinlP^l. In two dimensions, study on 
finite size system does not disclose any singularity at U c ^. In this paper, we focus 
on the single-site entanglement entropy S near the MT in the fermionic Hubbard 
model in infinite spatial dimensions. In this limit, the spatial fluctuations of elec- 
trons are suppressed while the local quantum fluctuation remains. The Hubbard 
model c an be solved exactly in this limit using the dynamical mean- fiel d th eory 
Recently, DMFT is used to evaluate the relative entropy^ as a 
measure of correlation for the Hubbard model as well as for a series of transition 
metal oxides^! For the half-filled Hubbard model, at low temperature, the Fermi- 
liquid state in small U regime is separa ted from the Mott insulator state in large 
U regime by a special second-order Qp r p2Z12212Q| Although this Mott transition in 
large spatial dimensions has received considerable attention in the past years, no 
analysis has been carried out for the critical behavior of its entanglement proper- 
ties. We expect that such analysis can help to deepen our understanding of the MT 
in general. 

The structure of this paper is as follows. In Section 2, we introduce the single- 
site entanglement entropy S and the method that we used to evaluate it for the 
Hubbard model, i.e., the DMFT with exact diagonalization and the two-site DMFT. 
In Section 3, we present numerical as well as analytical results for S near the MT. 
In Section 4 we end with a brief summary. 
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2. Model and Method 

The Hamiltonian of the Hubbard model reads (H = 1) 

i,j i ia 

Here, t%j is the hopping matrix element and U is the on-site repulsion of electrons 
with opposite spin. a a and c\ a are annihilation and creation operators of the elec- 
tron on site i with spin a, respectively, p, is the chemical potential. 

To study the local entanglement entropy, we divide the whole lattice into two 
parts, subsystem (A) (a single site i) and the environment (B) (the rest part of the 
lattice). For a given quantum state \tp) of the whole system, the reduced density 
matrix of the subsystem is 

p A = TtB|M# (2) 

The bipartite entanglement entropy between the subsystem and the environment 
is defined as (setting k = 1) 

S = -Tr(/5 A log 2 Pa) = -Tr(p B log 2 p B )- (3) 

From Eq.(2), one gets (i/j\Oa\iI>) = Tt(paOa) for any given operator Oa of the 
system. Taking Oi — 1, n^, n^, and n^Tiii for site i, one gets 4 equations about 
the diagonal elements of pi under the basis set (|0), | f}> | t-l})- Here we study 
the symmetry unbroken ground state of Hubbard model \ip). The off-diagonal ele- 
ments are all zero due to the U(l) and SU(2) symmetries of Hubbard model, i.e., 
tyWio-W) = (^14^4^) = (V ; |cl t c l 4,|'0) = 0. We therefore obtain 

/((l-n t )(l- % )) \ 

(n t (l-n ; )) 

((1-ntW) 

V (n tni )J 

The averages are with respect to the ground state of Hubbard model, and transla- 
tion invariance is assumed here. For a half-filled lattice in the paramagnetic phase 
where the Mott transition occurs, one has (n-j-) = (nj.) = 1/2. The entanglement 
entropy thus reads 

S = -2 Q - £>) log 2 Q - Dj - 2D log 2 D. (5) 

Here D = (n^ni) is the expectation value of the double occupancy. For the the 
ground state of the Hubbard model in infinite spatial dimensions, this quantity can 
be readily evaluated from the converged self-consistent solution of DMFT. 

The DMFT is a well-developed theory for treating the Hubbard-type strongly 
correlated models^! In DMFT, the Hubbard model is first mapped into an effective 
Anderson impurity model, 

H imp = Y [ e ka{ a a ka + T4 {a\ a c a + c^a^ + U 4 c t c l c ^ - ^ c ^ c °" ^ 

ka ' o 



PA 



(4) 



I 



I 
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Here a^a is the annihilation operator of bath site k, and the the parameters of the 
electron bath {ek,Vk} determine the dynamical "Weiss field" 

V 2 

G^iiuJn) = ioJn + /i - 22 ' ■ ( 7 ) 

lUJ n £fc 

k 

The impurity model is then solved to generate the impurity Green's function G(iuj n ) 
on the Matsubara frequency axis. Finally, through the self-consistent equation 

G{iu n ) = ^ . - Po{€) w , r de (8) 

J_ QC lUJ n +fJL — €— L(lUJ n ) 

together with Y>(iuj n ) = t/ ( 7 1 (ic<j„) — G~ l (iu n ), a new "Weiss field" Qq {ioj n ) can 
be obtained and used to update the bath parameters {e^, T4}. This process iterates 
and the converged solution of the impurity self-energy will be taken as the local 
self-energy of the lattice model, i.e., Ey(ia; n ) = Y^ imp (iui n )5ij . 
We adopt the semi-circular density of states 

2 



W — 1 is set as the energy unit. Eq.(9) is the density of states of free electrons on 
the Bethe lattice with infinite coordination number. It is widely used in the study of 
the MT because it simplifies the self-consistent equation while keeps the qualitative 
physics intact. We first solve DMFT equations using the exact diagonalization (ED) 
method of Caferral et a/PH Then, we resort to the two-site DMFT^l f or analytical 
results. At half filling, this theory is reduced to the linearized DMFtJSSI It is shown 
that it can produce rather accurate U c as well as physical quantities near the critical 
point. 

In the following, we present our results for D and S as a function of U, both 
from DMFT with ED and from the analytical two-site DMFT formulism. 



3. Results and Discussions 

3.1. Exact diagonalization results 

It has been know n ' | ou | that for finite temperatures < T < T c , the metallic state 
for small U is separated from the insulating-like state for large U by a finite regime 
U c i(T) < U < U C 2{T), in which the metal and the insulator phases coexist. The 
true first order phase transition occurs at U C {T) ( U C \{T) < U C (T) < U C 2(T) ) 
where the free energies of the metallic and of the insulating solutions coincide. At 
zero temperature, U cl w 2.38W^ and U c = U c2 « 2.94W^§1. The Mott transition 
becomes a special second-order phase transition at U C 2(T = 0): dE g /dU is contin- 
uous but the meta-stable solution of insulator extends from large U regime into the 
regime U cl {0) < U < U c2 (0). 

In Fig.l, we show the double occupancy D and entanglement entropy S as 
functions of U, obtained from DMFT with ED for a very low temperature T = 
lO -4 !^, which is practically same as zero temperature. As shown in Fig. 1(a), the 
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Fig. 1. (a) Double occupancy D and (b) entanglement entropy S as functions of U, obtained from 
ED calculation with N s = 6. Black solid dots are obtained by scanning from small U to large 
U , and red empty squares are obtained by scanning from large U to small U. Arrows mark the 
transition point U c 2- Other parameters are /x=U/2 and T = 10 _4 VF. 



double occupancy D decreases linearly as U increases up to the critical point U C 2- 
Near U C 2, coexistence of two solutions and hysteresis in {/-scanning are observed. 
From the solution of D obtained by scanning from small U to large U (solid circles 
in Fig. 1(a)), we can identify a linear form in the U < U C 2 regime after a non-singular 
term D c is subtracted, 

D - D c oc U c2 - U, (U<U c2 ). (10) 

The other solution of D (empty squares in Fig. 1(a)) is obtained by scanning from 
large U to small U. It extends to U c \ and recovers the first solution through a finite 
jump. From the continuation of D at U C 2 in the second solution, we can infer that 
in the U > U C 2 regime, a linear behavior with a much smaller slope must hold. 
Therefore, we can summarize D(U) as 

D = D c + a{U -U c2 ), (11) 

with a < and having different values on two sides of the phase transition. This is 

lo A I .... 

consistent with the scenario of the special second-order MT in infinite dimensions: 
d 2 E/dU 2 ~ dD/dU is discontinuous. Our ED calculation gives U c \ « 2 AW and 
U C 2 « 2.9W, in agreement with previous results of ED and projected self-consistent 
technique There is a small but finite double occupancy D c w 0.02 at U = U C 2- 



10, 2013 2:54 WSPC/INSTRUCTION FILE ws-mplb 



6 Dan-Dan Su, Xi Dai, Ning-Hua Tong 

This reflects that even at critical point and in the insulating state, there is residual 
local charge fluctuations. This has important consequence for the entanglement 
entropy in the insulating state. See below. 

In Fig. 1(b), the entanglement entropy S is shown. At U — 0, S — 2 comes from 
the equal population of electrons on the four local states. S decreases monotonously 
with U. After a inflection point at U = U C 2, it continues to decrease towards its 
strong-coupling limit Sqo = 1, which comes from the spin two- fold degeneracy in 
the paramagnetic insulator phase. Being consistent with D, two solutions coexist 
in the regime U c \ < U < U C 2 and a hysteresis is observed. What is interesting is the 
nonzero critical value of the entanglement entropy S c — Sao ~ 0.24. This reflects 
the residual entanglement between a single site and the others at the critical point 
and in the insulating phase. The critical behavior of S may be easily obtained from 



Here a < and it has different values on two sides of MT. Due to the finite D c , 
S is not singular at U = U C 2 but its derivative is discontinuous. This is the main 
result of this paper. 

Here is the big difference between the static mean-field theory and the DMFT. 
For static mean-field theories such as the Weiss mean-field theory, the entanglement 
will become zero as soon as the system enters the long-range ordered phase, while 
it stays as a constant S = 1 in the paramagnetic phase, arising from spin two-fold 
degeneracy. For DMFT, although the short range spin-correlation is not taken into 
account, the local temporal fluctuations do generate entanglement between one site 
and the others. Therefore, we conclude that for the Hubbard model, the local charge 
fluctuations in the insulating phase is associated with a nonzero local entanglement 
between one site and the others. 

Here we did not quest for the highest precision in our ED calculation. Near U C 2, 
the critical slowing down in solving DMFT equations prevents ED from obtaining 
rigorous conclusion. To get a more explicit solution, we resort to the two-site DMFT, 
which can give analytical results in the critical regime. 

3.2. Two-site DMFT result 

In the two-site DMFT, the lattice Hamiltonian is first mapped into an Anderson 
impurity model with one bath site, 



Eq.(8) and (5): 




(12) 




(13) 



The bath parameters e and V are determined by^l 



V 2 = zM$\ 



(14) 
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Here rii mp and ni at are electron density for the impurity site and for the lattice site, 
respectively, z is the quasi-particle weight and M2 is the second order moment of 
the free density of states, = pq(uj)uj 2 dtu . For the semicircular density of 

states in Eq.(9), = W 2 /A. For half-filling, the particle-hole symmetry in our 

model guarantees that e = 0. V is to be fixed by the second equation of Eq.(14). 

Eq.(13) and (14) can be solved numerically. The ni mp and ni a t are obtained 
from numerical diagonalization of Hi mp and from the lattice Green's function, re- 
spectively, z is obtained from the weight of the quasiparticle poles of the impurity 
Green's function. The numerical results for D and S are shown in Fig. 2. 

The two-site DMFT can also be solved analytically. For this purpose, we first 
solve the retarded Green's function of the impurity model at zero temperature 
fa = 0+), 

|<G|ct|n)| 2 IH4|G>I 2 



G a (u + ir,) = 



oj + if] + E n — E g oj + if] + Eg - E n 



(15) 



For the particle-hole symmetric case, the eigen values and eigen states of Hi mp can 
be solved analytically. The Green's function has four poles. Two of them near ±^[7 
are for Hubbard bands, and the other two poles near ±V 2 /U are precursors of the 
coherent Kondo resonance. To calculate z, we take the weights of the latter two 
poles and expand z in terms of powers ofV/U, 

36V 2 1584V 4 
z = -^--^ + .... (16) 

This expression together with Eq.(14) gives the solution for V as 

U 



V = —y/i-XP/Ul (U < U c ), (17) 



and V = for U > U c . Here U C = 6\J M 2 (0) . For the Bethe lattice, U c = 3W is very 
close the the ED results U C 2 ~ 2.9W. For U < U e , V > and the system is in the 
metallic phase. For U > U c , V = and the system is in the insulating phase. Note 
that this expression differs from that of PotthoflP^ in the prefactor, due to different 
ways of calculating J^S. 

Exact solution of Hi mp gives the analytical expression for the double occupancy 

as 



(U - VU 2 + 6W 2 ) , s 

D imv = n+ni = — i ; (18) 

P S T 11 4 (C/ 2 + MV 2 - UVU 2 + 64V 2 ) V ' 

Combining it with Eq.(17), one obtains the critical expression for D as, 

D imp = ^(l-U/U c ) (U<U C ), (19) 
and Di m p = for U > U c . For S we have 



s mp = ~ 1 1 - £\ fog 2 (i - ) (?■-.. r,.i. (20) 



i 



i 
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Fig. 2. Double occupancy D (a) and entanglement entropy S (b) as functions of U, obtained from 
numerical solution to the two-site DMFT equations. The solid lines are results for lattice model, 
and the dash-dot lines are for impurity model. 



and Simp = 1 for U > U c - 

However, it is noted that Si mp in Eq.(20) is not the correct results. Although the 
two-site DMFT equation requires rii mp = ni a t, the double occupancy Di mp does not 
equal to its lattice counterpart D lat , due to the fact that the two-site DMFT is an 
approximation to the full DMFT. In Fig. 2(a), we plot Di mp and D\ at as functions 
of U . Di a t is calculated using the lattice Green's functional It is seen that for 
U < U c , both D imp and D[ at are linear near U c . For U > U c , although D imp = 0, 
Di a t is small but finite, with a finite slope at U — U c , consistent with ED results. 

Correspondingly, S lrap calculated from Di mp becomes 1 immediately after U > 
U c , while Si a t calculated from Di at is not singular at U c . It has a linear behavior at 
U c but with a much smaller slope than that on the U < U c side. This agrees with 
the ED results. Therefore, we believe that the Si a t is an improvement over Si mp 
which ignores the finite D at U c . This is expected because Si a t takes into account 
the lattice information more elaborately than S, mp through the lattice Green's 
function. In conclusion, both the Si a t from two-site DMFT calculation and that of 
DMFT with ED give the critical behavior of the single-site entanglement entropy 
as described by Eq.(12). 

Our conclusion Eq.(12) is exact for the Hubbard model in infinite dimensions. 
Being different from the one-dimensionaff^ and two-dimensional casei^, our study 
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shows that MT in infinite dimensions fulfills the rigorous theorem by Larsson and 
Johannesson, which states that the discontinuity in the (fc — l)-th order deriva- 
tive of S gives a fc-th order QPT^, with k = 2 in this case. Besides the single-site 
entanglement entropy, two-site von Neumann entropy has been studied at the quan- 
tum phase transition^. For Hubbard model, it can be calculated using the cluster 
extension of DMFT^l 

4. Summary 

To conclude, we have analyzed the critical behavior of the single-site entanglement 
entropy S at the MT in infinite spatial dimensions, using DMFT with ED and the 
two-site DMFT. Even in the insulating phase, temporal local fluctuations generate 
nonzero S besides the contribution from spin two-fold degeneracy. The critical be- 
havior of S at the MT is not singular due to finite double occupancy at U c , but 
dS/dU is discontinuous at U c , being consistent with the second order scenario of 
MT and the Larsson- Johannesson theore rrP. 

This work is supported by National Program on Key Basic Research Project 
(973 Program) under grant 2012CB921704, and by the NSFC under grant number 
11074302. 
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